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Abstract — For a multi-user interference channel with multi- 
antenna transmitters and single-antenna receivers, by restricting 
each receiver to a single-user detector, computing the largest 
achievable rate region amounts to solving a family of non- 
convex optimization problems. Recognizing the intrinsic con- 
nection between the signal power at the intended receiver and 
the interference power at the unintended receiver, the original 
family of non-convex optimization problems is converted into a 
new family of convex optimization problems. It is shown that, 
for such interference channels with each receiver implementing 
single-user detection, transmitter beamforming can achieve all 
boundary points of the achievable rate region. 
Index terms — Gaussian interference channel, achievable rate 
region, beamforming 

I. Introduction 

The interference channel (IC) models a multi-user com- 
munication system in which each transmitter communicates 
to its intended receiver while generating interference to all 
unintended receivers. Determination of the capacity region of 
an IC remains an open problem. To date, the best achievable 
rate region was established by Han and Kobayashi in [1], 
herein termed the HK region, which combines rate splitting 
at transmitters, joint decoding at receivers, and time sharing 
among codebooks. The HK region was later simplified in [2]. 

Recently, [3] showed that the HK region is within 1-bit of 
the capacity region of the Gaussian IC. The results in [4]- 
[6], whose genie-aided approach is largely motivated by [3], 
established the sum-rate capacity of the Gaussian IC in noisy 
interference: treating interference as noise at both receivers is 
sum rate optimal, i.e., each receiver should simply implement 
single-user detection (SUD). In addition, even if the noisy 
interference condition is not satisfied, practical constraints 
often limit the receivers to implementing SUD. For example, 
the receivers may know only the channels associated with their 
own intended links. 

We assume in the present work that each receiver im- 
plements SUD, i.e., it treats interference as channel noise. 
In a preliminary work [7], we showed that beamforming is 
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optimal for the entire SUD rate region for a two-user real 
multiple-input single-output (MISO) IC. This result was used 
in [8] to characterize the beamforming vectors that achieve the 
boundary rate points on the SUD rate region. Later, the result 
in [7] was also used in [9] to derive the noisy-interference 
sum-rate capacity of the symmetric real MISO ICs. We note 
that the proof in [7] is applicable only to a two-user real MISO 
IC. 

There have been various studies concerning throughput 
optimization in a multi-user system under the assumption that 
each receiver treats interference as channel noise. However, 
even for the simple scalar Gaussian IC, computing the largest 
achievable rate region with SUD at each receiver is in general 
an open problem [10]. Exhaustive search over the transmitter 
powers is typically unavoidable due to the non-convexity of the 
problem. The difficulty is much more acute for the MISO IC 
case as one needs to exhaust all covariance matrices satisfying 
the power constraints, which renders the computation highly 
intractable. In this paper we propose an alternative way to 
derive optimal signaling for the SUD rate region for MISO 
ICs. Our approach is to convert a family of non-convex 
optimization problems for the original formulation to an equiv- 
alent family of convex optimization problems. What is more 
significant is that, given that each transmitter uses Gaussian 
input and each receiver implements SUD, all boundary points 
of the rate region can be achieved by transmitter beamforming. 

The rest of the paper is organized as follows. In Section II, 
we prove that beamforming is optimal for the SUD rate region 
of a multi-user MISO IC. Numerical examples are provided 
in Section III. We conclude in Section IV. 

Before proceeding, we introduce the following notation. 

> Bold fonts x and X denote vectors and matrices respec- 
tively. 

• (-) T and (•)T denote respectively the transpose and the 
Hermitian (conjugate transpose) of a matrix or a vector. 

• I is an identity matrix, is an all zero matrix, and 
diag(- • • ) is a diagonal matrix with its diagonal entries. 

• X > means that X is a symmetric positive semi- 
definite matrix. 

> tr(X) and rank(X) denote the trace and the rank of matrix 



X respectively. 

• (x)i denotes the ith entry of vector x, and X mx „ means 
that X is an m x n matrix. 

• ||x|| is the norm of a vector x, i.e., ||x|| = Vxhc. 

• E[] denotes the expectation. 

II. Multi-user MISO IC with single-user detector 

We define the received signal for user i of an m-user MISO 
IC as 



Yi 



T x- 



Ni. 



1, 



,m, 



(1) 



where ajj is the transmitted signal vector of user i with dimen- 
sion ti\ Yi is the scalar received signal; the Ni is unit variance 
Gaussian noise; and hji is the U x 1 channel vector from the 
jth transmitter to the ith receiver. The power constraint at 
the transmitter is tr(Sj) < Pi, where Sj = E[xixf]. We 
assume that receiver i knows only channel ha, and decodes 
its own signal by treating the interference from all other users 
as noise. The boundary points of the achievable rate region 
for this channel is characterized by the following family of 
optimization problems: 



i=l 



Hi Hi 



u- n 1, /, h^Sihii 

subject to IU = -log 1+ m " rr^T 



tr(Si) < Pi, Si h 0, 



1, 



TO, (2) 



where < /i, < oo. 

Apparently problem (O is a non-convex optimization prob- 
lem. For each choice of fi = [/ii,-- - ,/%], all possible 
[Si, • ■ • , S m ] must be exhausted. To obtain the entire SUD 
rate region, one has to go through this exhaustive search for 
all the fj, vectors. 

Following the same problem reformulation procedure used 
in [7], to characterize the achievable rate region of m-user 
MISO IC, it is equivalent to solve the following family of 
convex optimization problems: 

max hf l S l h ii 

subject to h, ,H,h,, < zfj, 

tr(Si) < Pi, S< hO 

i,j = l,...,m, i^j, (3) 

where zfj is a preselected constant denoting the interference 
power at the jth receiver caused by the ith transmitter. The 
problem reformulation is summarized in the following lemma. 

Lemma 1: For any vector fi with non-negative components, 
the optimal solution S* for problem (fJJ is also an optimal 



solution for problem <[3j with z. 



Proof: Problem (O is equivalent to the following optimiza- 
tion problem for the same fif 



max 




, < z* 2 



subject to h, S,h, 

ti(Si)<Pi, Si to, 



(4) 



The equivalence is due to the following. First, the maximum 
of problem (O is no smaller than that of problem (0]), since 
problem © has extra constraint hJjSihij < z* 2 . (This 
constraint is active, since the rates associated with S^ are on 
the boundary of the SUD rate region.) On the other hand, the 
maximum of problem (O is no greater than that of problem 
©, since the S*'s are also feasible for problem @. Therefore, 
problems (f2| and (0| are equivalent. We now recognize that 
problem (01 is equivalent to problem (f3]) by setting Zij = z*j . 

■ 

We remark that the optimization problem (0|i can not be 
solved directly as the constraint parameters z*^ all depend on 
the unknown optimal covariances. That is, unless the optimal 
S* of problem d2|i is obtained, the equivalent optimization 
problem © cannot be parameterized. Although we do not 
give an explicit solution of problem (f2]i for a given fi vector, 
Lemma[T]provides the following essential fact which is enough 
to obtain the entire SUD rate region defined by problem ©: 
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where the left-hand side denotes the collection of the optimal 
solutions of problem (fJI by exhausting fi, and the right-hand 
side denotes the collection of the optimal solutions of problem 
(O by exhausting zu. Since the SUD rate region is determined 
by the left-hand side of (O, Lemma 1 successfully converts a 
family of non-convex optimization problems (O into a family 
of convex optimization problems (|4j. 

Based on Lemma Q] we obtain the following theorem. 

Theorem 1: For an m-user MISO IC, the boundary points 
of the SUD rate region can be achieved by restricting each 
transmitter to implementing beamforming. 

Theorem [T] can be readily extended to complex channels. 
Before proving Theorem [T| we first introduce the following 
lemma. 

Lemma 2: Let x and y be two vectors with dimensions t\ 
and <2 respectively, and K >z be a (t± + £2) x (ti + t^) 
matrix with tr(K) < P. If 
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and Ku >_ is a preselected t\ x t\ matrix, then 

T 

K 



< 



x T K n x+ y/P-HiKn))' 



(6) 



(7) 



and the equality can be achieved by choosing K = K* 
1) When x T K n x ^ and \\y\\ ^ 0, we have 
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2) When i T Kni = and ||y|| 7^ 0, we have 

V^-tr(Ki- 



K* 



K, 



-K U W 



y/P-tv(K Ll ] 



P-tr(Kn) T 



where 1q 
with 
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being the eigenvalue decomposition and A being a strictly 
positive diagonal matrix. 
3) When |y|| = 0, we have 



K* = 



Kxi 




For all three cases, we have 

rank (K*) < max {rank (K n ) , 1} , 
The proof is omitted due to the space limitation. 



(10) 



(11) 



Lemma|2]is useful for the following optimization problems: 

T 
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, m, 
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max 
subject to 



,71, 

, m, 



subject to hi (Kn) = 0, i = 1, 
.</.. (Kn) < 0, 3 = 1, 
tr(K) <P, K y 0, 

where hi(-) and gj(-) are fixed functions. By Lemma |2] we 
can convert the above problem into 

( y/x T K u x + \\y\\ y/P - tr(K u 

/i i (K u )=0 ! i = l,. 

ft(Kn)<0, j = l, 

tr(Kn) < P, K b 0. 

Problems dT2b and ( fT3l ) have the same maximum. Once the 
optimal Kn for problem ( fT3l is obtained, we can construct 
the optimal K for problem ([T2]i by ([SJ, © and dTOj. We note 
that the choices of (O and (fT0T > are not unique. One can choose 
K different from that of © and (0 and still achieve the same 
maximum of problem (fT2b . 

With Lemma |2] we prove Theorem Q] as follows. 

Proof: By symmetry, it suffices to show that for the mth 
user, the optimal covariance matrix for the following 
optimization problem satisfies rank(Sj^) < 1: 

iT 



(13) 
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3 = 1, 
h 0, 
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where all the h m j's are £ 



x 1 vectors. 

We first show that problem (fl"4l i can be written as 



subject to hj Snhj < z m j, 



tr(Sn) 



< Prr 



j = 1,- • • ,m- 1 
Six ^0, 



(15) 



where /i and all the hj's, j = !,■■■ , m— 1, are m x 1 vectors, 
Su is an m x m matrix, and m is defined as 



min {t m , m — 1} . 



(16) 



Obviously, when m = i m < to — 1, problem ( TT4| ) is exactly 
problem (TT~5T > by choosing h = h mm , Su = S and hj = h m j- 
We need only to show the equivalence of problems (TBI and 
([T5l l when 772 == 771 — 1 < t m . 

Let the singular value decomposition (SVD) of h m i be 



\h rn l I 



(t m -l)Xl 



ftml — Ui 
and define S$ = Uf S m U a 

h mj= V l h mj, j = l,---,m. 
Substituting dT7]> and dT8j into (O, we obtain 
max fc (1)T SW/i (1) 

lnc " "'mm °m "•mm 

subject to ^B T SW^ < z 2 -, j = 2,- • 



(17) 
(18) 



,m— 1, 
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tr (SW) 



< p„, , siV y 0. 
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Consider h^ 2 an d let 
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elements of h^ 2 . The SVD of h m2 



is a vector consisting of the last t, 
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Define S£P = 
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On Substituting ([2Qj, (f23]) and (EU) into (O, we have 

2 max fc( 2 ) T s( 2 )/j (2) 



subject to fcS T S«)*« < 



j = 3,' • • ,777- 1, 



||ftml|| 
0(t m -l)xl 
(1) 
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S (2) 



< zt 



where Sn is an (to — 1) x (m — 1) matrix. The quadratic 
constraints in problem (|27T i are 
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= ftf Snft, < z 2 mp j = 1, • • • m - 1. (29) 



tr (s£>) < P m , S# b 0. 



(25) 



We note that the above transformation does not change the 
form of the existing constraint (see the third lines of problems 
(fT9l and d25ll). Now we continue the above procedure up to 
h m , m _i. In the jth transformation, we keep the first j — 1 
elements of ■ ^ and apply the SVD to the remaining (t m — 
j + 1) elements, and update the optimization problem. We 
formulate the jth iteration, j = 2, • • • , m — 1, as follows: 



Therefore, the quadratic constraints in problem d27T i are related 
only to Sn. By Lemma [2] problem (l27T i is equivalent to 
problem ([TBI . 

We summarize that we have shown the equivalence of 
problems ( fT~4T > and ([T5T l with all the vectors in (TT~5T > being mxl 
and Sn being to x to. 

By Lemma |2] we can reconstruct S in a way such that 
rank < maxjrank ^Snj , 1}. Let S^ be optimal for 
problem (1271 1. To prove Theorem [T] it is equivalent to prove 



rank 



< 1. 



(30) 



h U ~ 1] 



1,- ,3-1 
3 



M _ 



1(3) 



l 0-i)x0"-i) 




L (i-i)x(J-l) 



l (j-i)x(j-i)) 



U, 



,3-1 







(tm-i)xl 



Furthermore, it suffices to prove that the rank of the optimal 
covariance matrix for the following optimization problem is 
no greater than 1: 

max h T Sn/i 
subject to hJSuhj < z^j, j — L ■ • ■ , m — 1 

tr(s n )<P, SnhO, (31) 



where 



P 



*(s; x ) 



< Prr 



(32) 



jO'-i) 



r (i-i)x(i-i)) 
U, 



The equivalence is due to the fact that the optimal Sn for 
problem (15[ is also optimal for problem (|3~TT > and vice versa. 
Since Sn >: 0, we can define 



Sn — B B, 



(33) 



where (h^ 1] ) 
of 



denotes the jth to the i TO th elements 
and its SVD is 



where B is an m x m matrix. Then we can rewrite problem 
(ED as 



f^- 1)N ) 
\ mj J 



j,-,tr, 



(26) 



where UjUj = I(t m - J+ i)x(t m -j+i)- 

Finally, we convert problem (fl4l i into the following form: 
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z m j, 3 = l,---m- 1, 
S >r 0, (27) 



max ||Bft|| 
subject to IIB/ijH 2 < z 2 nj , 3 = 1, 
tr (B T B) < P. 

The Lagrangian of problem d34i > is 



L = -\\Bh\\ 2 +Y / ^{W Bh 3 

j=i 

+X m [tr (B T B) - P] . 

dL_ 
dB 



, m — 1, 
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Let 



B (C + A m I) = 0, 

m— 1 



where /i and h,- are (to — 1) X 1 vectors and h is a (t m — to + 

2 2 - 2 

1) x 1 vector. Furthermore, ||/i mm || = ||/i|| + /i . Let 



where C = —hh T + XjhjhJ 

3=1 

= H * diag [-1, Ai, • ■ • , X m -i] * H T 
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where H = [h,hi, ■ ■ ■ ,h m -i] is an to x to matrix, and C is 
an to x TO matrix. 

We then introduce the following lemma which is an exten- 
sion of Sylvester's Law of Inertia. 



Lemma 3: [11, Theorem 7] Let H be an m x n matrix 
and A be an n x n Hermitian matrix. Denote ir(-) and 
v(-) respectively as the numbers of positive and negative 
eigenvalues of a matrix. Then we have 

7r (HAH t ) < ir (A) , v (HAH t ) < v (A) . 

By Lemma [3] and the Karush-Kuhn-Tucker (KKT) conditions 
that require Aj > 0, i = 1, • • • , m, we have 



vr(C)<m-l, v(C) < 1. 



(38) 



Since C is an m x m matrix, we can write the eigenvalue 
decomposition of C as 



C = Q T diag(7/i,- • • ,?/™)Q 



(39) 



where Q T Q = I, and rji's are the eigenvalues of C in 
ascending order. From (|38l , i]i < and r)j > 0, j = 2, • • • , rn. 

Since the optimal B* satisfies tr (B* T B*) = P, from the 
KKT conditions we have A,„ > 0. Thus, we have 



rank (C + A m I) 

= rank (diag [A m + 771, X m + 772, • 
> m- 1. 



Since B is an m x m matrix, from 
optimal B* for problem (134-b satisfies 



' ' , A m + J/VnJ) 

(40) 

we conclude that the 
(41) 



rank(B*) < 1. 
Therefore 

rank (S ? *„) = rank (s^) = rank (B* T B*) < 1. (42) 

■ 

Remark: Theorem Q] proves the sufficiency of transmitter 
beamforming for achieveing the SUD rate region. However, it 
does not mean that the SUD rate region can only be achieved 
by beamforming. This depends on how we construct S after 
we obtain Sn. In the proof we choose S as ((H), © and ( TTOb in 
Lemma[2] However, the choices of (0 and ( TTOb are not unique. 
Another observation is that only (0 and (0 use full power. 
Equation ( TTOb corresponds to the case that h mm is linearly 
dependent of h m j, 3 ' — 2, • • • , m, so that h — in ( f27l i. This 
agrees with the result for two-user scalar Gaussian IC in which 
the maximum SUD sum rate sometimes is achieved when one 
user is silent [12, Theorem 6], 

III. Numerical examples 
Using Theorem [U we obtain in Fig. Q] the SUD rate region 
of a three-user MISO IC with the power constraint Pi = 1, 
P2 = 1.5 and P3 = 2. The channels are 
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where Hi = [hii,hu,hi 3 ], H 2 = [/i2i, ^22,^23] and H 3 = 
[/131, ^32, /133]. The solid curves are the rate regions for one 
user being inactive or at the maximum rate. That is, they are 
projections of the 3-D rate region onto a 2-D plane with one 
rate fixed at a constant value. 



IV. Conclusion 

We have considered multi-user MISO ICs where each 
receiver is limited to single-user detection. By exploiting the 
relation between the signal power at the intended receiver 
and the interference power at the unintended receiver, we 
have converted the original family of non-convex optimization 
problems into an equivalent family of convex optimization 
problems. Transmitter beamforming is shown to be sufficient 
to achieve all boundary points of the SUD rate region. 
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Fig. 1. The SUD rate regions of a three-user MISO IC. 
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